Optical and Thermal Transport Properties of an Inhomogeneous d-Wave 

Superconductor 
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We calculate transport properties of disordered 2D d-wave superconductors from solutions of 
the Bogoliubov-de Gennes equations, and show that weak localization effects give rise to a finite 
frequency peak in the optical conductivity similar to that observed in experiments on disordered 
cuprates. At low energies, order parameter inhomogeneities induce linear and quadratic temperature 
dependencies in microwave and thermal conductivities respectively, and appear to drive the system 
towards a quasiparticle insulating phase. 



PACS numbers: 74.25.Fy, 74.25. Gz,74.40.+k,74.80.-g 

Introduction. The study of disorder in high T c su- 
perconductors (HTSC) remains an engaging topic for 
at least two reasons: first, it is apparent that signif- 
icant disorder — perhaps originating with charge-donor 
impurities — is present in nearly all HTSC samples and, 
second, the controlled doping of substitutional impuri- 
ties is a powerful means of studying the electronic state. 
Transport experiments on the cuprates at low T have 
given us information on the quasiparticle lifetime in the 
disordered superconductor and indicated the existence 
of strong, near-unitarity limit scattering potentials as- 
sociated with impurities in the Cu02 plane. The sim- 
plest BCS-quasiparticle theories have been successful at 
describing qualitative features of transport experiments, 
but fail to explain many of the details. This is the prin- 
cipal motivation for the current work addressing trans- 
port properties of dirty d-wave superconductors. The 
approach taken here is to model the paired state as an 
inhomogeneous superfluid via the Bogoliubov-deGennes 
(BdG) equations. We focus most of our attention on opti- 
mally doped superconductors at low temperatures, where 
inelastic processes freeze out[Q and mean- field theory is 
most applicable. 

It is well known that HTSC are strongly affected by 
disorder because of the <i-wave symmetry of the pair or- 
der parameter Ay (i and j are site indices of the paired 
electrons). In a pure sample, the density of states (DOS) 
p{oj) is gapless and vanishes as \lu\ at the Fermi energy 
(taken to be here). A single strong-scattering impu- 
rity produces a pair of subgap resonances at ±wo (ojo < 
A mQX , A max the DOS peak associated with the gap edge 
in tunneling experiments). When a finite concentration 
m of impurities with impurity potential U is present, the 
isolated resonances are split, and broaden into an "impu- 
rity band" centered at the Fermi energy. The energy scale 
7 of the impurity band, below which p(E) crosses over 
from linear (\E\ > 7) to constant (\E\ < 7), is determined 
by m and U. These essential features are captured in 
the widely-used self-consistent T-matrix approximation 
(SCTMA) for impurity scattering §. The SCTMA is a 
perturbative scheme which is useful for treating point- like 
scatterers. It correctly incorporates physics associated 



with strong scattering potentials, but ignores correlation 
effects between impurities (i.e. localization effects), as 
well as the local response of the superfluid to the impuri- 
ties. In the SCTMA, 7 is also the quasiparticle scattering 
rate in the impurity band. 



While the SCTMA-based notion of an impurity band 
appears to be fairly consistent with the observed thermo- 
dynamic properties of the optimally doped cuprates ||, 
the simplest theory based on this picture disagree with 
transport experiments. Notably, the low-temperature be- 
havior of both the thermal and microwave conductivities 
in several systems disagree with the simple prediction, 
a, k/T ~ T 2 0. For example, the low-temperature mi- 
crowave conductivity in YBa 2 Cu307_,5 appears to vary 
roughly linearly with temperature, a ~ T^. In addi- 
tion, the optical conductivity of disordered cuprates is 
observed to have a maximum at a disorder-dependent 
frequency of order 100 cm -1 ||; this feature is also not 
found in the simple SCTMA analysis. @ Finally, the 
SCTMA predicts the universality of residual transport 
coefficients, i.e. limiting values of k/T and a as T — ► 
which depend only weakly on disorder S. While this has 
been confirmed in thermal conductivity measurements 
on YBa 2 Cu 3 07-,5§ and Bi 2 Sr 2 CaCu 2 8 [PL there are 
other systems where universality is not seen.|llj| We show 
below that some of these discrepancies can be understood 
within a BCS framework by going beyond the SCTMA. 



Approach. The BdG equations will be solved at two 
levels of approximation. Like the SCTMA, non-self- 
consistent (NSC) solutions assume that Ay is homoge- 
neous, but unlike the SCTMA, NSC solutions incorporate 
quantum coherence (ie. localization) effects associated 
with scattering from multiple impurities exactly. Self- 
consistent (SC) BdG solutions involve a further step in 
which the nonlinear response of Ay to the local disorder 
potential is determined. In both cases, the BdG equa- 
tions are solved on a tight-binding lattice with TV = 1600 
sites and up to 50 disorder configurations. In matrix 
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form, the mean-field Hamiltonian is 



* = 5>t 



At, 



(1) 



with <3>j = (c]|,Cij.). The subscripts i and j refer to site 
indices, and tij = —tSuj\ + (£7$ — fJ,)Sij with 5uj\ = 1 
for nearest neighbour sites, and otherwise. All ener- 
gies in this work are measured in units of t, and the 
lattice constant is a = 1. The bond order-parameter is 
Ajj = —V(cjiCij) with V the nearest neighbour pairing 
interaction. The pure ci-wave superconducting state oc- 
curs in the disorder-free limit and is related to the bond 
order parameters by Ay = ^A [(— l) Xij — (— l) Vij ] with 
( x ij ! Uij) connecting sites i and j, and where Ao is the ho- 
mogeneous d-wave amplitude. Spatial fluctuations arise 
naturally when one solves Ay self-consistently in the 
presence of a disorder potential. For this work V = 3.28, 
making Ao = 0.8 which is a factor of ~ 4 larger than the 
realistic case. The eigenstates are found using standard 
linear algebra routines to diagonalise Eq. 0. The quasi- 
particle DOS is p{u) = N- 1 {Y, n 5{uj - E n )) where E n 
are the eigenenergies for a given impurity configuration 
and (...) represents configuration averaging. 
The complex conductivity is 



a{w,T) = 



\ n,n 



i2 f( E n 

' flLO+ - 



■ f(E n >) 
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where 7" TO = (n\(p x /m) £§> T a \m) is the matrix element 
of the velocity between eigenstates n and m, r a (a = 
0, . . . ,3) is the Pauli matrix in particle-hole space, and 
oj + = to + i0 + . We use a binning procedure to evaluate 
the real part of a[w,T). Note the expression (|^) is not 
manifestly gauge invariant, but we do not expect this 
to be a problem since the collective response of the one- 
component charged order parameter occurs at the plasma 
frequency. 

For numerical reasons, we are restricted to eval- 
uating <j(oj,T) in the "gapless regime", 7 > T, 
which has been studied in intentionally damaged sam- 
ples of YBa 2 Cu307_(5^, [l2|]. In this regime, the 
SCTMA predicts that the real part of the conduc- 
tivity is <Jsctma{u -* 0, T) = <t o + a(T/~f) 2 and 
&SCTma(u,0) = (Too + o^'i^/l) 2 ! with a universal value 
for the residual conductivity: coo = j (j^Tiv a). In 
this expression, vp is the Fermi velocity and i>a = 
I V/c Afe I is the quasiparticle velocity component paral- 
lel to the Fermi surface. Both impurity vertex correc- 
tions and Fermi-liquid corrections renormalise (Too m an 
approximation |l3| where s-wave scattering is generalized 
to include anisotropic components, but where weak local- 
ization corrections and order parameter inhomogencitics 
are neglected. In this scheme the thermal conductivity is 
not renormalizedp^l , n(T)/T = HQQ + a(T/j) 2 , with uni- 
versal value koq — hk%[v F + vj^/vpv^ which survives 
this class of perturbative corrections. 
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FIG. 1: Finite frequency conductivity for fi = 1.2, T = 0, and 
A = 0.8 (ie. «f/«a = 2.5). (a) Curves are for self-consistently 
determined order parameter (SC BdG), homogeneous order- 
parameter (NSC BdG) and SCTMA with m = 0.04. (b) 
Density of states for NSC BdG and SCTMA with U = 10, 
and SC BdG with U = 5 (unitarity limit). The gap-edge 
for the tunneling density of states is A max ~ l.,the SCTMA 
impurity bandwidth is 7 fa 0.45. (c) Scaling of the NSC BdG 
conductivity for U = 10 and m = 0.02, 0.04, 0.06, 0.08, and 
0.14. Corresponding 7 are 7 = 0.29, 0.43, 0.55, 0.67, 0.91. A 
small deviation arises for m = 0.02 from finite size effects. 



Results. Typical results for p(u>) near unitarity — 
defined by c^o *C 7 and corresponding to U sw 10 and 
U ~ 5 in the NSC and SC calculations respectively — 
are shown in Fig. [I]. NSC calculations agree semiquan- 
titatively with SCTMA calculations except below an ex- 
ponentially small energy The SC result, by con- 
trast, shows a large disorder-induced suppression of the 
DOS relative to the SCTMA plateau. As discussed 
elsewhere Jl4| the "disorder-induced pseudogap" (DIP) at 
the Fermi-energy appears to be a generic feature of SC 
BdG solutions and has an energy scale related to luq but 
which grows with increasing m. We stress that for typ- 
ical planar Cu substituents in HTSC, this is an energy 
scale which is comparable to those explored in transport 
experiments. 

Figure shows the basic low-T result for the conduc- 
tivity with strong scattering impurities. <Jsctma{u) is 
approximately Drude-like for oj > 7, but saturates at 
the universal value (7], || (Too for u < 7. Numerical so- 
lutions of the BdG equations deviate significantly from 
this with o(uj < 7) linear in frequency, and rising to a 
peak at w « 7. This is true for both use an d &nsc 
(the SC and NSC BdG conductivities respectively), and 
is therefore the result of weak localization corrections 
to the SCTMA result Jl|. Indeed Fig. 0(c) shows that 
5<j(u>) = o nsc (^-0 ~ csctma{u) satisfies a scaling rela- 
tion, <5(t(cj)/<7oq = F{u>/^f), which is similar to the weak- 
localization scaling relation for dirty 2D metals. In con- 
trast, asci 1 ^) does not display a simple scaling relation, 
as we discuss below. At this point, we simply remark 
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that asc{u) always has less finite-w spectral weight than 
osctma^) for the same value of U [as illustrated in 
Fig. 0(a)], with the lost weight appearing in the super- 
fluid response [fl6|. 

The finite-frequency conductivity peak exhibited in 
Fig. is reminiscent of conductivity measurements in dis- 
ordered HTSCQ. Experimentally, the peak is also seen 
in the normal state T > T c , generally at higher energies, 
whereas in our approximation it occurs at significantly 
lower energies, and is much less pronounced than in the 
superconducting state. It is clear that inelastic scattering 
is important in the normal state and needs to be incor- 
porated in a complete explanation of the finite-frequency 
peak. Nevertheless, this is the simplest way of under- 
standing this feature of the optical data, which has not 
been reproduced in any other approach to our knowledge. 

The temperature dependence of the low- frequency con- 
ductivity is also of experimental interest. In Fig. we 
have plotted a(uJi,T), where uj\ = 0.0297 is the lowest 
nonzero frequency, chosen since uj = suffers from finite- 
size effects. The strong T-dependence of ctnsc(ui, T) 
is similar to the SCTMA result. On the other hand, 
asci^ijT) has a linear- T conductivity, reminiscent of 
most microwave conductivity experiments 0|. These re- 
sults are generic for a wide range of U near unitarity. 
Figure || also shows the thermal conductivity 



n(T) = 



1 



2-KhT 



dx : 



(3) 



where 



S T (x) 



2*2 



2n 2 h 
N 



n.n' 



\(v g ) n „'\ 2 6(x-E n )6(E n -E n ,), 



(4) 



with quasiparticle group velocity v g — + 



and 



(i/h)(xi — Xj)Aij in the site representation J 



For a finite system, the 5-functions in Eq. (|4|) are broad- 
ened to smooth the discreteness of the energy spectrum. 




FIG. 2: T-dependent conductivity, normalised to coo for (a) 
NSC BdG and (b) SC BdG with n, = 0.06, U = 10, fx = 1.2, 
and A = 0.8. U)i = 0.0297 is the lowest nonzero frequency 
used in calculating a with Eq. (^); k/(TLo) o(T) are evalu- 
ated using Eqs. ([|) and (|^). 



It is apparent in Fig. || that the Wiedemann-Franz law 
k/oT = Lq = fc^7r 2 /3e 2 , which is already violated be- 
cause of differences between the group and Fermi velocity, 
appears to also be violated at low T (at least in the NSC 
case) by weak localization corrections. For comparison, 
we show a similar calculation of the charge-conductivity 
which becomes exact in the limit u> <C T, 



dx 



9f(x) 
dx 



(5) 



with S a identical to St with the replacement of v g by 7 . 
From the figure, it is clear that Eq. (||) is in good agree- 
ment with a(u)i, T), and that both k(T)/T and cr(T) ex- 
hibit a linear T-dependence over a wide range of tempera- 
tures. The extent of the linear regime is discussed below. 
Linear power laws have been claimed in thermal con- 
ductivity measurements]^] (but other power laws have 
also been reported), and this work provides a potential 
mechanism. Conductivities with odd power-laws in T are 
difficult to achieve in the SCTMA because all quantities 
in the gapless regime are analytic functions of lu. We 
note that quasilinear behavior over some intermediate 
temperature regime has been obtained in SCTMA ap- 
proaches, by invoking a special combination of unitarity 
and weak scatterers[jl7|, by fine tuning scattering phase- 
shifts |l8|, or by including order parameter fluctuations 
in a SCTMA-like approximation fl9| . 

Finally, in Fig. || we study the dependence of crsc(u) 
on impurity concentration. As rii is increased, the peak 
position in use increases, in qualitative agreement with 
the scaling of <jnsc{w). The scaling of asci 1 ^) is not 
straightforward, however, since low-frequency spectral 
weight is depleted as increases, in contrast to (Jnsc{^) 
which depends only weakly on m at low w. The de- 
pletion is correlated with the growth of the DIP, shown 
in Fig. 0(b), reminiscent of disordered interacting met- 




FIG. 3: Scaling of asc{u) with rii at T — 0. (a) crsc(u) for 
a range of m between 0.02 and 0.14, (7 = 5. (b) Density of 
states at rii = 0.04 (open circles) and m = 0.14 (filled circles), 
(c) Thermal conductivity vs. temperature for same m, U. 
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als near the metal-insulator transition. For the model 
interaction chosen, we never observe a transition to a 
truly gapped state, and it is doubtful that a quasiparti- 
cle metal-insulator transition could be observed in real 
HTSC since superconductivity is destroyed in heavily 
damaged samples. However, the current work is strongly 
suggestive that the physics of disordered cuprate super- 
conductors is influenced by proximity to such a transi- 
tion. 

In Fig. ||(c), we plot the SC thermal conductivity 
for several impurity concentrations to illustrate the ro- 
bustness of the (quasi)linear-T regime. The regime is 
bounded by two disorder-dependent temperature scales; 
the upper crossover temperature is readily apparent up 
to n% — 0.08, and is correlated with the weak localiza- 
tion scale 7, while the lower bound signals a downturn in 
k/T which appears to scale with the DIP. For m — 0.14, 
there is no clear distinction between these scales. With 
current system sizes it is difficult to determine the lowest 
energy behavior. If one assumes the matrix elements of 
(v g ) have only weak energy dependence near the Fermi 
surface, then we expect St(x) ~ p{x) 2 which is ~ x 2a 
near x — (p ~ x Q [jl4j) in the SC BdG calculations. 
For sufficiently small T, then, one anticipates a down- 
turn with k/T ~ T 2a below the DIP energy scale. It is 
clear that a strong suppression relative to the universal 
SCTMA result k(T) — > koo is to be expected. 

Conclusions. We have observed effects with two dis- 
tinct physical origins in this work. First, there is a pro- 
nounced peak in a{u>) at u> ~ 7 arising from localization 
physics. This occurs whether or not the BdG equations 



are solved self-consistently, and is consistent with the ex- 
perimental fact that the peak is only observed in very 
disordered systems. Since 7 ~ y/nl for strong scatter- 
ers, our work suggests that a systematic study of sam- 
ples with varying impurity concentrations will provide 
an experimental means to distinguish between the weak 
localization mechanism presented here, and other pro- 
posed origins for the peak. Second, we have found that 
important physics associated with the correlated order- 
parameter response to disorder arises at low energies. 
Perhaps the most striking result is the observation of 
linear-T power laws in the charge and thermal conductiv- 
ities, as observed in some high-T c systems. In addition, 
order parameter supression effects appear to eliminate 
the residual conductivities at asymptotically low temper- 
atures expected on the basis of SCTMA and other treat- 
ments. We note that the current work assumes fairly 
disordered systems, and the extrapolation to the clean 
limit (7 < T) is not obvious. On the other hand, this is 
the first time that this additional source of off-diagonal 
scattering has been correctly accounted for in a transport 
theory, and there appears to be no reason in principle 
why these effects should also not be important in clean 
2D systems and possibly even in higher dimensions. To 
compare directly with experiments, the effect of realis- 
tic Dirac cone anisotropics and inelastic scattering needs 
to be better understood. Work along these lines is in 
progress. 
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